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1. Abstract : 

The quantum mechanical tunneling through multiple quantum barriers is a long-standing and well- 
known problem. Three methods proposed earlier to calculate the tunneling probabilities and energy 
splitting: (1). Instanton Method (2) WKb Approximation (3) Numerical Calculation.lnstaton 
method is helpful to understand the physical insight of quantum tunneling but the validity is 
restricted to the case of large separation between the two potential minima. WKB approximation is 
widely used in its simple mathematical form, but the result is inaccurate due to its inherent defect in 
connection formula. Recently WKB approximation has been developed by changing the phase lose 
at the classical turning points but no above approximation have provide the perfect result to the best 
of knowledge of Author. Using numerical methods, one can get the solution up to the desired 
accuracy, but a considerable deal of physical insight is lost in this process. In this paper, the Author 
presented the developmentof models of multiple quantum wells or barriers potential by using 
analytical Transfer matrix method (TMM), which has been applied to any arbitrary potential wells 
and barriers successfully. The author applied the above theory to three electronic device models and 
got satisfactory results. 

2. Introduction: 

Low dimensional carrier systems in the semiconductor heterostructures are gaining much 
importance in recent times due to the potential use of their unique properties in applications ranging 
from optoeltronics to high speed devices [1-4]. In this connection perpendicular transport of the 
carriers in semiconductor heterostructures has attracted much attention [5-8]. The MQW structures 
in particular, are becoming very important due to their potential use in the design and fabrication of 
quantum cascade lasers, resonant photo detectors, resonant tunneling diodes, single electron 
tunneling transistors [8] etc. Moreover, with the decrease in the dimensions of the CMOS devices 
the effect of tunneling of carriers becomes very important in order for estimating the various 
leakage currents flowing through the devices present in the VLSI chips. 

In this paper an attempt has been made to study the tunneling of carriers through a quantum barrier 
and plot the variation of the transport coefficient with respect to carrier energy. The range of energy 
include the classically forbidden transitions as well. The transmission/tunnelling coefficient, which 
is the flux of particles penetrating through the potential barriers to the flux of particles incident on it 
at the other interface, has been computed by using the Ben Daniel-Duke (BDD) boundary 
conditions [9] for solving the Schrodinger equations for the electrons (carriers in this case) inside 
the coupled well regions and the barrier in between. The theory has been based on the transfer 
matrix method. Tunneling depends significantly on the barrier width. Scaling of structure 
dimension affects this variation very sharply. 
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The material pairs of interest include CdS/CdSe, AlGaAs/GaAs and InAs/AlSb. InAs/AlSbs based 
HEMTs are excellent for application in satellites due to their low operating voltages [3]. Devices 
based on AlGaAs/GaAs have been in use over quite some time and CdS/CdSe QW structures 
promise of improved gain performance for light emitting applications [2] 

The effective masses of the carriers are different inside the well and in the barrier, which are made 
of different materials. Further, this effective mass may change with energy as is given in the case of 
the InAs/AlSb pair [10]. This affects the tunneling behaviour of the electrons because there is a 
dependence of the transmission coefficient of the carrier effective masses. In this paper the effect of 
variation in effective mass on transport coefficient has also been studied. The variation of transport 
coefficient with electron energy has been studied and compared for different barrier widths for each 
of these material pairs. 

Low dimensional carrier systems in the semiconductor heterostructures are gaining much 
importance in recent times due to the potential use of their unique properties in applications ranging 
from optoeltronics to high speed devices [1-4]. In this conection perpendicular transport of the 
carriers in semiconductor heterostructures has attracted much attention [5-8]. The MQW structures 
in particular, are becoming very important due to their potential use in the design and fabrication of 
quantum cascade lasers, Quantum Infra-red photo-detectors (QWIP), Quantum cascade detectot 
(QCD), resonant photodetectors, resonant tunneling diodes, single electron tunneling transistors [8] 
etc. Recently much research is carried on in the field of semiconductor infrared detectors. These 
are based on inter-subband transitions and have the possibility of device tailoring for distinct, 
narrow band detection energy range. Quantum cascade detectors (QCD)s are based on a bound-to- 
bound transition. Here the excited state is coupled to an extraction cascade / phonon-stair that 
transport electrons vertically to the ground state of the next period that follows. [9] 

In this paper an attempt has been made to derive a general mathematical model that will help in 
evaluating the energy Eigen values inside asymmetric and aperiodic MQW structures. The model is 
tested on the experimental data for a 16.5 fim QCD of reference 9. This model will also help in 
studying the tunneling of carriers one well to another through a quantum barrier. The transmission 
coefficient, which is the flux of particles penetrating the potential barrier to the flux of particles 
incident on it at the other interface, has been computed by using the Ben Daniel-Duke (BDD) 
boundary condition [9] for solving the Schrodinger equations for the electrons (carriers in this case) 
inside the coupled well regions and the barrier in between. The theory has been based on the 
transfer matrix method [10-16]. Tunneling depends significantly on the barrier width. Scaling of 
structure dimension affects this variation very sharply. 
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The structure on which the mathematical model is tested is an array of In 5 3 Gs.4 7 As/rn 52 A1 48 As 
layers of thicknesses in Angstrom in the order 157/55/69/26/90/43/94/40/95/38/96/28 [17]. The first 
two wells have widely different widths and the barriers that separate them and flank the pair on 
either side are of unequal width. The effective mass of the carriers are different inside the well and 
in the barrier, which are made of different materials. The model is tested on this well pair to match 
the computed result of coupled energy states with those obtained experimentally. The agreement is 
satisfactory. 

Inside the structure the wells and barriers of widths 90/43/94/40/95/38/96 can be approximated by a 
periodic pattern as the wells widths are varying very slightly and so do the barrier width. Under this 
circumstance the general model can be simplified and the energies are computed. The comparison 
with experimental data prove encouraging. 



The quantum mechanical theory of tunneling through a classically forbidden energy state can be 
extended for other types of classically forbidden transitions. In this paper a generalized theory of 
quantum tunneling for transition through multiple quantum barriers has been developed and 
tested for three different pairs of materials. These adjacent lower energy regions, that are 
separated by a quantum barrier, are coupled and this is the general pattern for many of the 
heterostructures. The electron wave functions in the lower energy regions and the barrier region 
are obtained by solving the Schrodinger equations that satisfy appropriate boundary conditions 
[11,12]. The solutions depend on the effective masses of the carriers in the regions concerned. 
Hence the tunneling probability will not only be a function of the dimensions of the barrier alone 
but will be affected by change of materials as well as its energy dependent effective mass. 



3.(a) Methodology - 1 Phase: 




Fig. 1. Schematic illustration of a multiple quantum wells (MQW) heterostructure. 
The total wave function 0(r,x) for the electron can be written as 
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<D(r, x) = exp (ik r) Y(x) 1 
Where T(x) satisfies the one dimensional Schrodinger equation with position-dependent electron 
effective Mass, x represents the growth direction of heterostructure. 

The Schrodinger equation according to the effective mass theory for the finite potential barrier 
and the well regions on either side take on the well-known form 



h 2 d 2x ¥(x) 



+ (E- V)^(x) = 



2m* dx 2 

with appropriate effective masses m* and potential energies V for the region where the equation 
is defined. Inside the barrier region the effective mass m*= m B and potential energy is Vo while 
those outside the barrier are respectively given by m*= m w and zero. 

Now, we use the transfer matrix for the jth junction and w should generalize the result for In- 
junctions. This matrix relates the coefficients of the wave function at one end of the junction to 
those of the other one, so that the wave function can be written as 

¥j w (x) =Aj W exp(ik w x) +Bj w exp(-ik w x) 3 
for the jth well, and 



¥j b (x) =Aj b exp(ik b x) +Bj b exp(-ik b x) 
for the jth barrier 



Where k v 



h 2 



and kj,= 



2m h (V -E) 

h 2 



with m w and m B are the electron effective 



masses in the regions outside the barrier and barrier region respectively. 

By matching the continuity of the wave function T(x) and its appropriate normalised derivative 
1 d y/(x) 

— at the boundaries and form 2x2 transfer matrix equations for each interface and we 

m dx 

derive a matrix formula that relates the coefficients Aj and Bj with 
Aj+i andBj+i. 



B 



M 21 [J] 



M 22 [J] 



v'+i . 



By using equation (5) we obtained the coefficients of the wave function at the leftmost slab to 
those of the right most slabs 
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Where M . is the jth transfer matrix corresponding to the jth junction written as: 

M, =M b {b J )M w (a j )M h (b j+l ) 7 
Where bj and are the widths of the jth barrier and jth well respectively, M 6 (6.)and 
M w (a . ) Correspond to the transfer matrices for the jth barrier and jth well respectively. 

Here the total transfer matrix is expressed as the cascading of a series of individual barrier and 
well. 

From equation (6) & (7) it is readily found the transmission amplitude Q is given by 

2 



M n + M 22 +i(k w M 12 -k w l M 2l ) 



Where M tj are the elements of the total transfer matrix. 



From the definition of transmission coefficient (T), we obtained the following expression from 
equation (8) 

T=|ef 9 



From equation (5) we obtained the expansion coefficient of the wave function at the left most slab 
to those of the right most slab as: 



= M [1] M [2] M [3] M [hl] 



A, 



=M 



total 







Which satisfies M 2 i = 
We obtained the energy Eigen values and Eigen functions by solving equation (12). 
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11 

12 



The wave vector k w outside the barrier is a real quantity for all positive energies E of the electron. 
The conduction band minimum in the region outside the barrier is taken as the zero of energy. 
Since the conduction band minimum of the barrier is above that of the region outside hence for 
certain energies of the electron (E<Vo) the wave vector k b will be imaginary and for energies 
above Vo (E>Vo) it will be real. Thus the energy is divided into two regions, (E<Vo) for non- 
classical transition by tunneling and (E>Vo) where transition can take place even under classical 
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conditions. The solutions will be different and the methods applied to evaluate the transmission 
coefficients in the two regions will be different. Since k w and k b are both dependent on effective 
mass and energy, hence for different material pairs the variation of the transmission coefficient 
for energy values normalized with respect to the barrier height will be different. [13, 14] 
Here P = m B / m w is called the mass discontinuity factor. It plays a very important role in 
determining the transmission coefficient. 

For perfect transmission through a barrier sandwiched between two wells we must have 
p (k b /k w ) - k w /(pk b )] 2 sinh 2 k w L = 13 
where L is the length of the barrier. So there will be transmission at energy values E n = Vo + 
(n7ih/L) 2 /2m B where transmission coefficient becomes 1 . The barrier becomes transparent at these 
energy values. This is called Ramsauer-Townse [13] effect in atomic physics. The resonance 
condition is satisfied at normalized energy values (E/Vo) = l+n 2 /p, where n is an integer. The 
resonance values therefore depend upon the mass discontinuity factor. 

When the electron energy E is less than the barrier height Vo the situation is somewhat different. 
The above derivations assume that the mass discontinuity factor remains unchanged with energy. 
However, it has been pointed out [15] that for InAs/AlSb material system the effective mass 
varies with the energy according to the expression 

m(E) = m*[l+(E-E c )/E eff ] 14 

With E c represting the conduction band minimum and E eff is the effective bandgap at the energy E. 
The transmission coefficient normally has a dependence on energy through the expression of the 
wave vector k b . The energy variation of the mass discontinuity factor P(E) adds another 
dimension to this variation. It is worthwhile to examine the effect of this variation and how far it 
affects the tunneling phenomenon. 

The transmission coefficient of electrons through a potential barriers is important for studying the 
leakage current in MOSFETs with dimensions in the nanometer range. It is also a crucial 
parameter for studying the behavior of multiple quantum well structures where the barrier is 
sandwiched between two coupled quantum wells. In the second case the equations are modified 
to include the well dimensions also. When both the well and barrier regions are in the nanometer 
range we expect further quantization of the energy levels. This is being considered in a further 
study of the multiple quantum well structure. 



3. (b) Methodology - 2 na Phase: 
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The quantum mechanical theory of tunneling through a classically forbidden energy state can 
be extended for other types of classically forbidden transitions. In this paper a generalised 
theory of quantum tunneling for transition through a barrier between two quantum wells has 
been developed where the electron wave functions in the well region and the barrier region are 
obtained by solving the Schrodinger equations that satisfy appropriate boundary conditions 
[10,1 1]. The solutions depend on the effective masses of the carriers in the regions concerned 
and the height of the barrier on either side of the well. Hence the tunneling probability will not 
only be a function of the dimensions of the barrier alone but will be affected by change of 
materials as well. 

We have started with the most generalised case where the well widths and barrier widths are 
all unequal. In addition to this the well and barrier materials can be different with different 
barrier hights and effective electron masses. If we assume that the barrier widths are such that 
adjacent wells are coupled through the barriers then we will have to deal with five wave 
functions; the wave functions of the electrons inside the wells and the electron wave functions 
for the barriers on either side of these wells. We consider the Schrodinger wave equations for a 
series of three barriers of widths 2d, 2b and 2f with wells of width 2a and 2c between the 
barriers 2d and 2b and 2b and 2f respectively. The well barrier combination is shown in figure 
1 below. 



2d 




2b 




2f 




2c 




2a 















-b-2e-2d -b-2c -b +b b+2a b+2a+2f 



Figure 1: Araay of MQWs. 

The Schrodinger equation according to the effective mass theory for the finite potential barrier 
and the well regions on either side take on the well known form 

(h 2 /2m*) d/dx {dWdx} +(E-V)T =0 1 . 

with appropriate effective masses m* and potential energies V for the region where the 
equation is defined. Inside the barrier region the effective mass m*= m B and potential energy 
is V while those outside the barrier are respectively given by m*= m w and zero. 
Solution of the Schrodinger equation yields the electron wave functions inside the barrier and 
outside which is represented as for the most general case: 

Tbi(x) = A B i exp (ik B ix) + B B i exp (-ik B ix) for; -(b+2c-2d)<x< -(b+2c) in the first barrier 

2. 
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1*wi( x ) = A W1 exp (ikwix) + B wl exp(-ik wl x) for -(b+2c) <x<-b 3. 
Inside the first well 

^B2(x) = A B2 exp(ik B2 x) + B 2 exp(-ik B2 x) for -b <x< +b 

in the second barrier 4. 

1*wi( x ) = A W 2 exp (ikw 2 x) + D 2 exp(-ik w2 x) for b<x< b+2a 

inside the second well 5. 

^ m x) = A B3 exp (ik B3 x) + B B3 exp (-ik B3 x) for (b+2a)< x<(b+2a+2f) 

in the second barrier 6. 

Where k Bn = (2m Bn [Vo n - E]/h 2 ) 1/2 and k W i = (2m W iE/h 2 ) 1/2 with m wi and m Bn are the electron 
effective masses in the well and barrier regions respectively. In the most general case the 
barrier heights Vo n and electron effective masses m wi and m Bn in the different regions are 
assumed to be all different. 

The boundary conditions that are applied assume continuity of the wave functions at the well- 
barrier boundaries i.e., T Bn (x) = Tw^x) at each boundary. 

The other important boundary condition called Ben Daniel Duke (BDD) boundary condition 
relates the flux of electrons incident on the barrier to that which penetrates it. This is expressed 
mathematically as 

(l/rn wi )[d¥ wi (x)/dx] = (l/m B „) [dT Bn (x)/dx] 

at well-barrier interface. 7. 

By matching the continuity of the wave function T(x) and its appropriate normalized 
derivative as described above the transfer matrix method can be applied to find out the 
transmission coefficient x = I(F/A) 2 I at the farthest region. The Matrix equation is very 
complex and involves a large number of terms. The coupling energy between two adjacent 
wells can be obtained by putting x = 1 and varying the energy value to obtain the best fit. This 
iterative method will yield the energy for which electrons may be coupled to both the wells by 
penetrating through the intervening barrier. 

This mathematical model is tested on an array of In 53 Gs 47 As/In 52 A1 48 As layers of thicknesses 
in Angstrom in the order 157/55/69/26/90/43/94/40/95/38/96/28 [17]. The first two wells have 
widely different widths and the barriers that separate them and flank the pair on either side are 
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of unequal width. The effective mass of the carriers are the same in both the wells 
(m w i=m w2 =m w ) but different from that in the barriers, which are made of the same material so 
that m B i=m B2 =m B3 =m B . The barrier hights in all the barrier regions are same and equal to V . 
The transmission coefficient derived by Transfer Matrix method for the most general case 
simplifies to some extent and becomes 

x =[ 4 k B k w /( m w m B )] 4 /[T 1 +2T 2 +2 T 3 +2 T 4 ] 

8. 

where 

Ti={(k B /m B ) +(k w /m w )} 8 + {(k B /m B ) -(k w /m w )} 8 + 6{(k B /m B ) 2 -(k w /m w ) 2 } 4 
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T 2 ={(k B / m B ) +(k w /m w )} 4 {(k B / m B ) 2 - (k w /m w ) 2 } 2 [cos(4{bk B +ck w })+cos(4{ak B +bk w }- 
2cos(4ck w ) - cos(4ak B ) - cos(4bk B ) - cos(4{ak w +bk B +ck w })] 

10. 

T 3 = {(k B / m B ) 2 -(k w /m w ) 2 } 4 [4cos(4ak w )cos(4ck w ) -2cos(4bk B ) { l+cos(4[a-c]k w )} -2 {cos(4ak w 
)+cos(4ck w ) } 1 -cos(4bk B ) } 11. 

T 4 ={(k B / m B ) -(k w /m w )} 4 {(k B / m B ) 2 -(k w /m w ) 2 } 2 [cos(4{bk B -ak w }) +cos(4{bk B -ck w }) - 
cos(4ak w ) - cos(4bk B ) - cos(4 {bk B -ak w -ck w })] 12. 

In order to compute the coupling energies the transmission coefficient is equated to 1 and 
energy values are obtained by iterative method. 

When the well widths and barrier widths are respectively same and the same material is used 
for the well regions and the barrier regions respectively then the equation for the transmission 
coefficient simplifies to equation 13. 

(m w k B )/(m B k w ) = tan(ck B )/tan(ck w ) where c is the well width. 

13. 

The energy values are computed by iterative method. The transmission coefficient takes the 
form given below. 

x = 16(k B k w ) 2 /[( m B m w ) 2 {(k B /m B ) +(k w /m w )} 4 + {(kw/m w )-(k B /m B )} 4 - 2{(k w /m w ) 2 -(k B / 
m B ) 2 } 2 cos(4ck B )}] 14. 

At the coupling energies obtained from equation 13 the transmission coefficients given by 
equation 14 will tend to unity. 

4. (a) Results: 1 st Phase: 

In this section we present our results obtained numerically by using MATLAB programming for 

the transmission coefficient across quantum wire containing multi -barrier heterostructure. The 

pairs of materials chosen are CdS/CdSe, AlGaAs/GaAs and InAs/AlSb. The parameters used in 
the computation are given in table 1. The effective mass of electrons in Al x Gai_ x As depends on 
the mole fraction of x, where x represents the concentration of Al. 

In case of A^Ga^As the effective mass is given by m A1GaAs = (0. 063+0. 083x) and the energy 
band gap is E A1GaAs = (1.9+0. 125x+0.143x 2 ). 

For each pair three values of the barrier width L are taken; 5nm, 1 Onm and 20nm. In all cases the 
transmission coefficient increases with diminishing dimension, as expected. This increase is the 
slowest for InAs/AlSb, which has the highest value of p\ The effective mass variation for 
AlSb/lnAs pair is included in the calculations. It is found that the value of the mass discontinuity 
factor does not change very sharply with energy for AlSb/lnAs. However this slow variation of 
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mass discontinuity factor P appears to slow down the increase of transmission coefficient, 
especially at very low barrier width (5nm). 
Table- 1: 



Parameter 


CdS/CdSe 


Al x Gai. x As/GaAs (x= 0.47) 


AlSb/lnAs 


Conduction 

Band-gap 

AE C 


Egcds: 2.36eV; 

Egcdse: 1.69eV; AE C = 67% of 

(Egcds-Egcd S e) = 0.45 eV 


EgAiGaA S : 1.99eV; 

EgGaAs- 1.42eV AE c =67%x 

(EgAlGaAs - EgGaAs) = 0.38 eV 


Ecoin As : 0.0 eV; 
EcoAisb: 2.1eV; 

AE C = (EcOAlSb - EcOinAs) 

= 2.1 eV 


Effective 
mass (m*) 


m*cds: 0.20m o ; 
m* C ds e : 0.1 3mo 


m* G aA S : 0.067m o ; 
m*AiGaA S :0.106m o 


m*i nAS : 0.020m o ; 
m^Aisb: 0.098m o 


Mass 

discontinuity 


p = 1.54 


P =1.58 


p = 4.9 



Quantum tunneling for all barrier widths is least for AlSb barrier in InAs/AlSb pair and most 
significant for AlGaAs barrier for the in GaAs/AlGaAs/GaAs. This appears to be quite justified 
because InAs/AlSb/lnAs structure has the highest barrier height and GaAs/AlGaAs/GaAs the 
lowest.(Figures 2, 3, 4) 

As the barrier width decreases the tunneling increases and transmission coefficient value rises 
with normalised electron energy. The rate of rise is sharpest for GaAs/AlGaAs/GaAs and rather 
slow for InAs/AlSb/lnAs. This is easily explained if the barrier heights are compared. 
CdS/CdSe/CdS structure lies midway. 

For (E nor = E/Vo) <1, the transmission coefficient increases from to 1 in a non-linear fashion. 
For each pair, the transmission coefficient is lower for wider wells as expected. However, the rise 
is the slowest for InAs/AlSb/lnAs because of the highest value of the mass discontinuity factor. 
Beyond the normalized energy (E nor = E/Vo) >1, there is resonance; i.e., there are quantised 
energy values where transmission reaches peak values sharply. This variation is prominent for all 

three structures, especially near E nor = 1. As the barrier width decreases the peaks get more 
separated in energy; this is most prominent in GaAs/AlGaAs/GaAs structure. For 
AlSb/lnAs/AlSb the difference between the maxima and minima of the values of the transmission 
coefficients remain relatively unchanged with decreasing well width. For CdSe/CdS/CdSe and 

GaAs/AlGaAs/GaAs the maxima of transmission coefficient gradually increase as E nor increases 
beyond l.The difference between the maxima and minima also gradually decreases. 
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Here the regions of lower band gap are taken to be semi-infinite. Hence the normalized energy 
values vary continuously. If the width of the lower band gap materials outside the barrier are 
reduced to the order of nanometers then the energy values inside the quantum wells will be 
quantised. This effect will be reflected in the nature of variation of the transmission coefficients 
with normalised energy and is expected to change significantly. This will have a crucial effect on 
carrier tunneling in multiple quantum well structures and leakage current in field effect devices. 
Studies of this effect are being explored by the author. 



CdSe/CdS 



L=5nm 



GaAs/AIGaAS L=5nm 



InAs/AISb 



L=5nm 





E/vO 



(a) 



E/vO 



(b) 



E/vo 



(c) 



Fig. 2: Variation of transmission coefficient of electrons with normalized energy E/Vo for (a) CdSe/Cds (b) 
GaAs/AlGaAs and (c) InAs/AISb for well with. 



CdSe/CdS 



L=10nm 




E/VO 



(a) 



GaAs/AIGaAs L=10nm 



hs.5 




E/vO 



InAs/AISb 



L=10nm 




E/VO 



(b) 



(c) 



Fig. 3: Variation of transmission coefficient of electrons with normalized energy E/Vo for (a) CdSe/Cds (b) 
GaAs/AIGaAs (c) InAs/AISb for well with 10 nm. 
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CdSe/Cds 



L=20nm 




E/vO 



GaAs/AIGaAs 
L=20nm 




0.5 1 1.5 



E/vO 



InAs/AIS 



L=20nm 




E/vo 



(a) (b) ( c) 

Fig. 4: Variation of transmission coefficient of electrons with normalized energy E/Vo for (a) CdSe/Cds 
(b) GaAs/AIGaAs (c) InAs/AlSb for well with 20 nm 



4. (b)Results: 2 nd Phase: 



The model is tested on a MQW structure obtained from a paper by F.R. Giorgetta et al [17]. 
The effective electron mass of the well material In.53Gs.47As is taken to be m w = 0.043mo [18] 
and that of the barrier layer of In 52 A1 4g As is m B = 0.072m o . The well and barrier widths that 
alternate are given in Angstroms as 157/55/69/26/90/43/94/40/95/38/96/28 [17]. The 
conduction band offset for In.53Gs.47As/In.52Al.4gAs pair is 0.51eV.[17] 

The array is divided in to two parts for testing the model. The first part consists of the first pair 
of wells of different widths lying between the three unequal barriers. The energies for which 
the transmission probabilities are unity are E^O.0205 eV, E 2 =0.1014eV, E 3 =0.1959eV, 
E 4 =0.2899eV and E 5 =0.3910eV. These computed values agree quite well with the 
experimental values given in figure 1 of reference 17. The experimental values are around 
0.02eV, O.lOeV, 0.20 eV, 0.30 eV, 0.40eV and 0.45eV. The values cannot be determined 
with greater accuracy as they are not quoted but shown as the possible Eigen values 
corresponding to the Eigen functions. Tunnelling at around 0.2eV is very distinctly indicated. 
In our computation the transmission coefficient is 0.99 at 0.205 eV. For the other energies 
transmission coefficients are high but less than what we have for 0.205eV. The transmission 
coefficients at these energies are around 0.97. 

The model is applied to the last three wells. Here the array of In.53Gs.47As/In.52Al.4gAs pair is 
90/43/94/40/95/38 Angstroms. The well widths of 90, 94 and 95 Angstroms are quite close. 
So are the barrier widths of 43, 40 and 38 Angstroms. Here the simpler expressions 
represented by equations 13 and 14 may be applied without much error. The energies 
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Computed are 0.78 eV and 0.205 eV. The transmission coefficients for tunneling approach 
0.99. 

5. Conclusions: 

In conclusion we may add that the general model developed can be used for designing MQW 
structures. To achieve tunneling through the structure at any specific energy value the well 
widths can be computed iteratively for x = lat the specified energy. Applying this model to 
design MQW structures will be done in future research work. 
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